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We propose the implementation of a strongly-driven one-atom laser, based on the off-resonant 
interaction of a three-level atom in A-configuration with a single cavity mode and three laser fields. 
We show that the system can be described equivalently by a two-level atom resonantly coupled to the 
cavity and driven by a strong effective coherent field. The effective dynamics can be solved exactly, 
including a thermal field bath, allowing an analytical description of field statistics and entanglement 
properties. We also show the possible generation of Schrodinger cat states for the whole atom-field 
system and for the field alone after atomic measurement. We propose a way to monitor the system 
decoherence by measuring atomic population. Finally, we confirm the validity of our model through 
numerical solutions. 

PACS numbers: 42.50.Fx, 42.50.Vk, 73. 21. La 



I. INTRODUCTION 

In cavity quantum electrodynamics (CQED) the inter- 
action between atoms and photons can be investigated 
experimentally under carefully controlled conditions, and 
described by relatively simple models [1[ . These features 
make CQED an almost ideal framework to investigate the 
foundations of quantum mechanics and their application 
to quantum information 0. For instance, two-atom en- 
tanglement |S| as well as the entanglement between an 
atom and a photon [3| have been recently demonstrated. 
On the other hand, the basic interaction between a two- 
level atom and a cavity field mode, as described by the 
Jaynes-Cummings (JC) model [11, leads to nonclassical 
effects carefully tested in recent years @. Furthermore, 
it allowed the implementation of the micromaser Q and 
the microlaser [Sj in the strong coupling regime of CQED, 
in the microwave and optical domain, respectively. Fur- 
ther efforts led to the implementation of a trapped ion 
as a nanoscopic probe of cavity field modes More 
recently, a single trapped neutral atom in a high-Q opti- 
cal cavity (lo| allowed the implementation of a one-atom 
laser i.e., lasing with only one intra-cavity atom. 

These systems can exhibit features that are not present in 
standard macroscopic lasers such as thresholdless genera- 
tion and sub-poissonian photon number distribution [l^ . 

Another milestone in CQED experiments was reached 
in Ref. , where a "Schrodinger cat" state of the cavity 
field, a mesoscopic superposition of two coherent states, 
was realized. There, the field decoherence was moni- 
tored through atom-atom correlation measurements p^ . 
State reconstruction of nonclassical intra-cavity fields 
was also possible through atom-cavity dispersive interac- 
tions [l^, [l^ . More recently, a remarkable proposal for 
the resonant generation of Schrodinger cat states was 
implemented in the lab [l8| , and tested with the help of a 
quantum spin-echo technique . The understanding of 



entanglement in atom-cavity systems was enhanced when 
an additional driving field acting on the cavity mode was 
added on top of the atom-cavity JC interaction [20l. |2]|. 
In this respect, recently, an elegant analysis of a driven 
cavity containing a two-level atom explained the absence 
or increase of ent ang lement in the transient of the atom- 
cavity dynamics [2^ . Unfortunately, most realistic mod- 
els including dissipative processes require numerical anal- 
ysis, or ideal theoretical conditions for the sake of semi- 
analytical derivations. 

In this work, we introduce an integrable model of a 
strongly-driven one-atom laser (SDOAL) operating in the 
optical regime of CQED, where the coherent driving field 
acts directly on the atom. We consider a realistic model 
consisting of a three-level atom in A-configuration placed 
inside a single-mode optical cavity, coupled off-resonantly 
to three coherent laser fields. We show that this model 
can be reduced to two atomic levels coupled to a cavity 
mode and a strong classical driving on the atom. In this 
strong-driving limit [23] , we solve the full system dynam- 
ics [24], in the transient and in the steady state, provid- 
ing one more of the few examples of an exactly solvable 
open quantum system. In previous works (25l . [26j . we 
developed related results for microwave cavity fields and 
two-level Rydberg atoms, not a good model for a field 
in the optical regime and fast decaying atomic dipolar 
transitions. Here, we solve analytically the master equa- 
tion for the full atom-field system in the SDOAL model. 
Next, we exploit the obtained solutions for the analysis 
of atom-field entanglement and the decoherence of atom- 
field superposition states ( Schrodinger cat states) via 
the measurement of atomic populations. In addition, the 
generation and the decoherence of cat states of the cavity 
field alone is described. Finally, we present numerical re- 
sults confirming the validity of the approximations made 
to derive the effective master equation of the SDOAL 
model. 
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The paper is organized as follows. In section II we in- 
troduce the integrable model of a SDOAL. In section III 
we solve analytically the master equation for the atom- 
cavity dynamics. In section IV we consider the dynamics 
of the cavity field and atom subsystems. In section V we 
describe entanglement properties and the environment- 
induced decoherence of the SDOAL, presenting a scheme 
to monitor decoherence via atomic populations measure- 
ments. A numerical analysis that confirms the validity 
of the model is presented in section VI. Conclusions are 
reported in section VII. 



II. THE STRONGLY-DRIVEN ONE-ATOM 
LASER MODEL 

We consider a three-level atom (ion) in a A- 
configuration trapped inside an optical cavity (Fig. [T]). 
We assume that the transition |1) |2) is quadrupolar 
and, hence, the metastable states |1) and |2) cannot be 
coupled directly, but only via the level |3). The level |3) 
can decay via spontaneous emission and, therefore, the 
external lasers and the cavity field are all far detuned 
with respect to the corresponding transition frequencies. 
We suppose that the atom interacts off-resonantly with 
a single mode of a cavity field of a frequency ujf on the 
transition |3) |2). The same transition is coupled 
off-resonantly to a coherent field of a frequency The 
remaining atomic transition |3) |1) is coupled off reso- 
nance to two lasers of frequency uj[ and w. The different 
frequency detunings, A and A' < A, of these two A- 
processes prevent the system from undesired transitions. 




FIG. 1: Atomic energy levels and the applied fields. A and 
A' denote the frequency detunings, a;/ is the frequency of the 
cavity mode, and lj'i, uj'2, uj are the frequencies of the lasers 
applied to the associated transitions. 

We assume, without loss of generality, that both 
the cavity mode coupling frequency g and the associated 
laser Rabi frequencies Q,,U,'i,U,'2 are real. The Hamil- 
tonian 'H{t) for the whole system can be written as 
ii[t) = 'Ho+Hi{t), where 

Tio = hujsS^^ + nw2S^^ + hLOiS^^ + hiufo^a, (1) 



+ Ml'^^e-^'^'^^Sf + e"^'--^S'^_^). (2) 

Here a {a^) is the cavity mode annihilation (creation) 
operator and, following the notation of [l^l, we define 
the atomic operators as follows, 

Sf = \3){2\,S^J = \2){3\,Sl' = \3){l\, 

S]^' = |1)(3|,5'" = |J)(J|, (J = 1,2,3). (3) 

We rewrite the Hamiltonian Ti.{t) in the interaction pic- 
ture leaving the unavoidable time dependence in the term 
related to the laser Rabi frequency fl 

+ hg[ia+^)Sf + h.c.] 
9 

+ ftr!'i[(l + ^e'(^-^')*)^f ^H.c.]. (4) 

If {A, A', |A - A'|} > {g, n, the unitary dy- 

namics of the atom-field system in Eq. (|4]) can be de- 
scribed by an effective Hamiltonian for a two-level atom 
coupled to the cavity mode and in presence of a classical 
driving field. This is due to the fact that, under these 
conditions, the energy diagram of Fig. [1] can be under- 
stood as composed by two independent A-schemes. In 
this case, it is straightforward to prove that we can build 
the second-order Hamiltonian 

n'^s = -hg.sia^Sl^ + aS'J) - ^17cff(^+ + S'J), (5) 

with gcs — gfl/A and = ^[^'2/ A'. In Eq. (O, as 
is usually done, we have assumed the compensation of 
constant AC Stark-shift terms by a proper retuning of 
the laser frequencies. The Stark-shift term depending on 
the intracavity photon number can be neglected if $7 ^ 
In the strong-driving limit, QcS 3> gcS, as explained in 
Ref. [2^, we can derive the final effective Hamiltonian 

As = -h~^{a^ +a){Sl^ + S'J). (6) 

We have tested numerically the above analytical consid- 
erations and proved, in fact, that Eq. ([6]) describes the 
correct effective dynamics. However, we want to show 
here that we can go beyond the limit of uncoupled A- 
schemes and obtain a similar dynamics with less demand- 
ing conditions on the experimental parameters. To prove 
this statement we exploit the small rotations method of 
Ref. ,23], which is essentially a perturbative method for 
deriving effective Hamiltonians. First, we introduce the 
operators of a SU (3) deformed algebra, 

123 ^ ^ ^)^23 ^ ^23 ^ (gt + ^^)523^ 

9 9 
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where bt = I + fp-e^^^ ^ Using the identity relation 

/ = S^^ + S^^ + S^^ we can rewrite the interaction Hamil- 
tonian Eq. ((4]) in the compact form, 

Hiit) = -HA' + hA'S^^ ~ hiA~ A')a^a 

+ hg{Xf + x^^) + hn[{Yl^ + Yl^). (8) 

We can ehminate the dependence of Tii on the upper 
level 1 3) by applying two consecutive small rotations. The 
first unitary transformation U^^ — exp{a{Y^^ — Yl^)}, 

with a = ^ <C 1 and the condition (^)^ ^ 1, al- 
lows us to eliminate the dependence on operators Y^'^. 
The second unitary small rotation, given by U^^ — 
exp{P{Xf - Xl^)}, with /3 = ^ < 1 and the condi- 
tions (^)^ ^ 1, (^^-)^ ^ 1, can be used to eliminate 
the dependence on X^?. After some lengthy algebra we 
derive the effective two-level Hamiltonian 

Ti'eft = h{A~ A')S^^ + h{A~ A')a''a 

- hg,s{a^Sl^ +aS'J)-hn,s{Sl^ + S'J), (9) 

where we introduced the effective coupling and driving 

frequencies geS = ^ and flcs = ^^r- Note that the 
above derivation does not depend on the order of the two 
small rotations. The effective Hamiltonian ([9]) is exact to 
zero-th order in the diagonal terms and to first-order in 
the other ones. From now on we shall consider the case 
of small detuning difference <C 1, such that the 

diagonal terms are negligible. This will allow us to ob- 
tain an exactly solvable model of system dynamics even 
in the presence of dissipation, that is a one-atom laser. 
Hence, the initial model described by the Hamiltonian 
of Eqs. Il]) and ([2]) reduces to a Hamiltonian that ex- 
hibits an effective coupling of the states |1) and |2) to 
the cavity mode in the presence of a classical external 
field driving the atomic transition. Now we can apply 
the unitary transformation U = exp{—i^leg{S^ + S]?)t} 
to obtain [23| 

H'W = _!_)(_!+ e-2'"-*|+)(-| 

- e2^""«*|->(+l)a^+H.c.], (10) 

where we used the eigenstates |±) = ^^^^^^ of the oper- 
ator Sx = + S]^. In this way, we put in evidence fast 
rotating terms in Eq. (fTO|) and, after applying the RWA 
with fioff 3> 5cff, we obtain the final effective Hamilto- 
nian 

= -n^ia^ +a)iSl' + S',^). (11) 

This Hamiltonian has the structure of resonant 
and simultaneous Jaynes-Cummings and 'anti-Jaynes- 
Cummings' interactions [l^ and its dynamics is better 
understood in terms of Schrodinger cat states than Rabi 
oscillations, as will be discussed later. Note that the 



Hamiltonian of Eq. ((TT|) is similar to the one of Eq. ([6]) 
but with a more relaxed set of parameters. Furthermore, 
the dynamics is fully confirmed by numerical simulations. 

To describe the open atom-cavity system dynamics we 
must include the dissipative effects due to the coupling 
of the cavity to the environment. We note that the de- 
cay of the upper level |3) can be neglected because of 
the elimination procedure described above. Therefore, 
the system dynamics can be described by the following 
SDOAL master equation (ME) 

1 ^ 

Paf = --^[HcB, paf] + CpAF, (12) 

where the dissipative term is the standard Liouville su- 
peroperator for a damped harmonic oscillator 

^PAF = -^{a^apAF - '^apAFa^ + PAFO^a). (13) 

Here, k is the cavity photon decay rate and we consider 
the limit of zero temperature because the system operates 
in the optical regime. 

III. ANALYTICAL SOLUTION OF THE SDOAL 
MASTER EQUATION 

The time evolution of the atom-field system is de- 
scribed by the density operator pAF^t) which is the so- 
lution of the ME in Eq. (fT^ . In order to solve it, we 
introduce the following decomposition 

PAF{t) = \+){ + \® Plpit) + \-){-\® p2F{t) 

+ \+){-\® P:iF{t) + \-){+\® pMt) (14) 

Here, piF{t) {i — 1,...,4 ) are operators describing the 
cavity field defined as 

pMt) = { + \pAFit)\+), p2Fit) = {-\pAF{t)\-), 
P3F{t) = { + \PAF{t)\-), p4F{t)^ {-\pAF{t)\+). 

(15) 

Then, the master equation (|12p is equivalent to the fol- 
lowing set of equations for the operators piF(t) 

/9l,2F = ±i^[a^ +a,pi^2F]+ Cpi^2F, (16) 

P3AF = ±'i^{a' + a,P3AF} + ^PSAF, (17) 

where brackets [ , ] and braces { , } denote the standard 
commutator and anti-commutator symbols. In order to 
describe a one-atom laser dynamics we assume that the 
initial atom-field density operator is paf{0) = ® 
|0)(0|. Therefore, operators piF{0) read 

P.F(0) = i|0)(0| (*-l,...,4). (18) 

This choice is suitable in the optical regime of CQED also 
because the generation of coherent states is difficult due 
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to the very fast decay of the cavity mode. Nevertheless, 
from a theoretical point of view and for an extension to 
the microwave regime of CQED it is possible to general- 
ize the following analysis to the case of a field prepared 
in a coherent state. The results in Eqs. (|3ip should be 
modified by redefining the form of the function a(t). 
In order to solve Eqs. (|16p and (|17p . we map them onto a 
set of first order partial differential equations for the func- 
tions x^{^^,t) = TrF[p^F{t)D{(3)], I = 1, ...,4, where i)(/3) 
denotes a displacement operator. The functions Xi(/3ji) 
cannot be interpreted as characteristic functions for the 
cavity field, because the operators pipit) do not exhibit 
all required properties of a density operator. As a conse- 
quence the functions XiiPi^) do not fulfill all conditions 
for quantum characteristic functions. Nevertheless, they 
are continuous and square-integrable, which is enough 
for our purposes. From Eqs. p6|) and (|17p we obtain the 
following set of partial differential equations. 



dx 



1,2 



dt 



(19) 



dx. 



3,4 



dt 



f(4+^*|.)X3,4. (20) 



To solve these differential equations we use the method 
of characteristics, for which it is useful to rewrite them 
in terms of the real and imaginary parts of the complex 
variable P — x + iy, 



dxi 



dt 



T«5cffa;Xi,2 

2(-'+2/^)Xi,2-2(-a^+ya^)Xi,2,(21) 



K , d d . 



(22) 



If in the equations for ^3,4 we introduce the shifted vari- 
able y — y T 2-^ the above equations can be written 
as 

%f + f(^^+^^)^^'^ = ^i,2(-,y)Xi,2, (23) 

where Hi^2{x,y) = x[F[{x) T ^^2(2;)] + vG'{y) and 
if3,4(x, y) = y[E[{y) T E'M] + xD'{x). There, we have 



also introduced the derivatives of the following functions 
Fi{x) = -■^a;^, F2{x) =ig^iix, G{y)^-jy'^, 
Ei{y) = --,y\ i?2(y) = 25,ff(^lny + y), 

4 K 



D{x) 



(25) 



With these definitions, together with the initial func- 
tions Xi,o(a;,?/) — Xii^jVi^) associated with the ones in 
Eq. (flS]) . we can write the time-dependent solutions as 



Xi,2(x, y,t)^l exp{-^^ T 2*^[1 - e-^]}, 

Z Z K 



(26) 



X3,4(a;,y,t) = — exp{ =F2-^[l-e 2]}, 



(27) 



where 



fit) = exp (-2^t + 4^(1 - e-"*/2)| (28) 
The most striking feature of the solutions for Xza in 



Eq. (p7)) is the presence of the factor e in contrast 

to the solutions for xi.2 in Eq. ((26|) . This factor leads to 
the vanishing of functions X3.4 for sufficiently long times. 
To better understand the solutions (pS)) and (P?]) we 
rewrite them in terms of the complex variable [3 

Xi,2(/3,t) = iexp|-^±[/3«*(t)-/3*a(i)]|, 

(29) 

X3.4(/3,t) = Z|lexp|-^T[/?a*(0+/3*«W]}, 

(30) 

where we have introduced the complex time-dependent 
function a{t) = i^(l — e"''*/^). We immediately recog- 
nize that the operators pip (t) corresponding to the func- 
tions XiW, t) are 



PMt) = -\ait)){a{t)\, 



P2F{t) 



a{t)){-a{t)l 



P^Fit) = l-II^\a{t)){-ait)\, 

P4F(i) = ^-^|-«(^))(«(t)|. (31) 

We describe now the generation of Schrodinger cat 
states for the whole atom-field system. Actually, in the 
limit of Kt <C 1, when the unitary dynamics dominates 
over the incoherent cavity dissipation, f{t) ~ e"^'"'*-" 
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so that the state of the atom-field system is well approx- 
imated by: 

\^{t))AF = ■^{\+)\a{t)) + \-)\- am, (32) 
with a(t) = i^. 

On the other hand, the steady state of the atom-field 
system is the mixed state 

paf - i+>(+iprF + i->(-iP2i^ 

|+)(+||a-)(a-| + |-)(-||-a-)(-«-| 



(33) 



with a^'^ = igcs/K. 



IV. CAVITY FIELD AND ATOM SUBSYSTEM 
DYNAMICS 

We consider the reduced density operator for the cavity 
field pFit) = Tr A[pAF{t)] = piFit) + p2F{t), where Tta 
denotes the partial trace over the atomic variables. From 
Eq. ([31]) we obtain 



PFit) 



y{t)){a{t)\ 



v{t)){-a{t)\ 



(34) 



and we see that it is always a mixed state. The cavity 
field mean photon number after an interaction time t is 



{N){t)=TrF[a^apF{t)] = \a{t)\' 



(35) 

In the steady state the cavity field mean photon number 
is given by (N)^'^ = (gcs/n)^, that is, the squared ra- 
tio between the effective coupling frequency and the cav- 
ity decay rate, which rule the coherent and incoherent 
regimes of cavity field dynamics, respectively. The time- 
dependent photon number distribution Pn{t) is given by 
a Poissonian distribution 



Pn{t) 



Ht)\ 



2n 



-i"(t)r 



(36) 



corresponding photonstatistics are 
pL'''Ht) - {n\p^F-'\t)\n) 



1 



|a(t)P 



1 ± f_i)» -^(^^ 



(38) 



In the transient dynamics, for times kt <^ I, a(t) = a{t) 
and the cavity field states are even and odd cat states 



\m)''F^^ 



\a{t))±\-a{t)) 

^2(1 ±e-2|"(t)P) 



(39) 



As is well known [28|, states as in Eq. can exhibit 
quantum effects including oscillating photonstatistics, 
sub Poissonian photonstatistics, and quadrature squeez- 
ing. In Fig. [5] we show the time behaviour of the Mandel- 

Fano parameter Q = ^ — — I in both cases of atom 

detected in the lower {Q^^\t)) and upper 

(Q(2)(t)) state 

and for different values of the steady state mean photon 
number. We see that Q(^'(t) and Q^^H*) exhibit super 
and sub Poissonian photonstatistics, respectively, before 
approaching the steady state Poissonian distribution. 



1 



Q 
0.5 



-0.5 



-1 
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FIG. 2: The Mandel-Fano parameter Q for the cavity field 
versus dimensionless interaction time in the case of atom 
detected in the lower (Q'^'(^)) ^"^^ upper {Q^'^\t)) states. 
We consider different values of the mean steady-state photon 
number {N)ss- (a) 1, (b) 2, (c) 5, (d) 10. 



Hence, at any time the photon number distribution of the 
SDOAL is that of a coherent field, a natural consequence 
of tracing orthogonal atomic states |±). Certainly, this 
will not be the case if we make a projective atomic mea- 
surement in the bare basis {|1), |2)} at a time t during 
the transient. Actually, after the atom measurement, the 
cavity field is in either of the pure states 



pr\t) 



plFjt) + P2F{t) ± [psFjt) + P4F{t)] 



(37) 



where Pi,2(0 the probability to find the atom in the 
state |1), 1 2) respectively at a time t (see below). The 



Now, we consider the reduced density operator for the 
atom pA{t) = TrF[pAF{t)]- From Eq. (pij) we derive the 
following density matrix in the basis {|+), |— )} 



P^it) = l 



1 fit) 
fit) 1 



(40) 



From this atomic density matrix we can derive the prob- 
abilities pi^2 to find the atom in the lower or upper state. 



1, 



pi,2(i) = (l,2|p^(t)|l,2) = -[l±/(t)] 



(41) 
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We observe that in the steady state the atomic popula- 
tion of the upper level |3) is zero and those of the lower 
and intermediate levels are both equal to 0.5. The phys- 
ical intuition behind this result is the orthogonality of 
coherent states \ct{t)) and | —a{t)) when t oo. We will 
employ these results in the following section to study the 
entanglement properties and the decoherence of the sys- 
tem. 



0.3r 
0.2^ 



0.1 




V. ENTANGLEMENT AND DECOHERENCE 
ANALYSIS 

In sections III, we presented a new scheme for gen- 
erating atom-field superposition states [see Eq. in 
the transient regime and we described the steady state 
of a SDOAL. Now, we evaluate atom-field entanglement 
properties and show how to monitor the decoherence to- 
wards a steady state. We have shown that the state of 
the whole atom-field system is almost a pure state on a 
time scale much shorter than the cavity decay time 1/k. 
Therefore, in this case, we can use the entropy of en- 
tanglement E{^) as an entanglement measure. It can be 
calculated in a straightforward way using the equality , 



Ei^) = Sa^Sf, 



(42) 



where Sa and Sp denote the von Neumann entropy of the 
atomic and field subsystems, respectively. The atomic 
entropy reads 



(43) 



where {Ai, A2} are the eigenvalues of the reduced atomic 
density matrix pAii)- In the limit nt <^ 1, the atomic 
density matrix in Eq. (|40p can be approximated by 



PA{t) = \ 



1 

-2|a(t)|2 



g-2|a(t)|^ 
1 



(44) 



whose eigenvalues are 



Ai,2(i) 



lie 



-2|Q(t)|- 



1 r 
2 



lie" 



V2] 



(45) 

In Fig. [3] we plot the time evolution of the von Neumann 
entropy Sa for different values of the steady-state mean 
photon number {N)^^ . We see that the system gets more 
entangled for larger values of {N)^^ , i.e., when the ratio 
gcs/k is large. 

The next question is how to monitor the decoherence of 
the whole atom-field system, that is, the reduction from 
a pure state to a statistical mixture. The environment- 
induced decoherence of a cavity field prepared in a su- 
perposition state has been both theoretically and exper- 
imentally studied in the case of high-Q microwave cavi- 
ties [1^, [l3|. We remark that the cavity field reduced 
density operator does not depend on the decoherence 
function f(t). However, a simple way to monitor the 




FIG. 3: The Von Neumann entropy Sa for the atom-field 
system versus dimensionless time for different values of the 
steady-state mean photon number {N)ss'- (a) 0.25, (b) 1, (c) 
5, (d) 10 , (e) 20. 

decoherence of the atom-field system is to measure the 
atomic populations pi, 2 (i) ofEq. (|41]) . In fact, the atomic 
inversion I{t) — pi{t) — P2{t) is exactly the function f{t) 
(Eq. (l28l) '). that can be rewritten as: 

f{t) = cyi^[-2{Nf^Kt + 4(7V)^^(1 - e-'=*/2)]_ (4g) 

In Fig. m we illustrate the time evolution of the atomic 
inversion I{t) for different values of (N)^^ showing that 
the decoherence dynamics is rather complex. 
For dimensionless times Kt 1 the inversion I{t) shows a 
Gaussian fall-off as exp{—{N)^^K?t^/2}, independent of 
the cavity field decay rate k. We recall that in this limit 
the interaction generates the atom-field cat-like superpo- 
sition as in Eq. ([5^ . After this transient, the effective 
decoherence process begins in correspondence to the in- 
flection point at time: 



2 

^-ln(l 

K 



1 + 16(N) 



ss 



8{N) 



ss 



(47) 



For (N)^^ ^ 1, corresponding to effective strong cou- 
pling conditions, we have tp ^ — . — gZ^ , and the 

decoherence function f(t) can be well approximated by: 
litp) e^p{-2K{Nf^il - e-^*^/^)it - tp)}. (48) 
Hence, we can introduce the decoherence rate 



7d 



1 + 16(iV)'5S - 1 



ss — 



Soft- 



(49) 



We note that the decoherence rate is given by the effective 
coupling constant and hence it is independent of cavity 
dissipation. 

In the opposite limit of small (N)^^ we have tp ^ 
k~^; the system is close to the mixed steady state and 
we recover an exponential decay with a decoherence rate 
7^ = 2k{N)^^ . This is the standard decoherence rate 
for cat states of the cavity field alone [l3| , which can be 
generated in our system after an atomic measurement. 
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function by means of the following main rule 



4 Kt 5 



FIG. 4: Atomic population inversion / versus dimensionless 
time evaluated for different values of {N)ss'- (a) 0.25, (b) 0.5, 
(c) 1, (d) 5, (e) 20. 

VI. NUMERICAL SIMULATIONS 



In this section we discuss the theoretical approach pre- 
sented in the above sections from a numerical point of 
view. In fact, by means of a first order perturbative ap- 
proach we have reduced the full three-level system dy- 
namics to an effective two-level one described by the ME 
in Eq. We now discuss the validity of that approx- 

imation for both Hamiltonian and dissipative dynamics. 

In the numerical analysis we need to solve the full sys- 
tem ME 



PAF = --^[Htit),pAF] 



Cp 



AF, 



(50) 



where the system Hamiltonian is given in Eq. ([5]) and the 
dissipative process is ruled by the Liouville super opera- 
tor in Eg. p^ . To solve Eq. ([50]) numerically we consider 
the dimensionless time t = At and the following dimen- 
sionless parameters 



A' 



g ~ n 

A' A' 



A' 



A ^ 



A 

(51) 



In order to solve the ME by means of the Monte Carlo 
Wave Function approach (MCWF) 29], we rewrite the 
ME in the Lindblad form to identify the collapse and the 
"free evolution" operators 



PAF = --{HePAF- PafHD + CpafC^ (52) 



where the non-Hermitian effective Hamiltonian Tig is 
given by 



(53) 



and the only one collapse operator is C* = \/%a. The 
system dynamics can be simulated by a suitable num- 
ber of trajectories, i.e. stochastic evolutions of the wave 



m+6i))^{ 



if Sp{i) > Nrnd 



(54) 

where dt is a suitable small time interval, 6p{t) is the col- 
lapse probability at time t, and Nmd is a random number 
generated from a uniform distribution in [0, 1]. We note 
that the collapse probability depends on the cavity field 
mean photon number {N){t) and can be evaluated as 
dp{t) — Stk{N){t). In the simulations we must consider 
parameters values in agreement with the theoretical con- 
ditions required by the two small rotations. 

First we discuss the full three-level system Hamilto- 
nian dynamics (k = 0) in order to confirm the validity 
of the effective two- level Hamiltonian of Eq. pT|) . We 
consider the time evolution of the cavity field mean pho- 
ton number and of the atomic populations, and we com- 
pare the numerical results with the theoretical expres- 

sions {N)it) = and pi^^it) = i[l ± exp(-^)], 

Psit) — 0. As an example, we show in Fig.[5]a case where 
the ratio of the effective parameters is ~ 25. We see 
a good agreement for the mean photon number (Fig.[S^). 
The theoretical functions pi,2{t) fit the envelopes of the 
numerical fast oscillating populations (Fig. [S)d). In fact, 
in the numerical analysis we do not take into account the 
RWA approximation. In particular, the populations of 
levels |1) and |2) approach the expected value of 0.5, and 
the population of the upper level |3) is always negligible. 
In addition, we tested the prediction that the effective 
dynamics allows to generate cavity field cat states when 
the atom is measured in level |1) or |2). In Fig. [6] we 
show the Wigner function that describes in phase space 
the cavity field state prepared by an atomic measurement 
in level |1) and we see the typical features of a cat state. 

Now we consider the full dynamics including dissipa- 
tion of Eq. ([50]) . for the same parameters as in Figs. [5l [SI 
and with k = ^eff , so that we expect that the steady state 
value of mean photon number is one and it is reached 
in a time that is twice that of the atomic populations. 
In FiglTK we compare the numerical results for the time 
evolution of the cavity field mean photon number to the 
theoretical behavior predicted by Eq. p5|) . showing that 
there is a good agreement. In Fig. [7)3 we consider the 
numerically simulated time evolution of the atomic pop- 
ulations Pj{t) {j = 1,2,3) compared to the theoretical 
functions in Eq. (|4T|) . We remark that the population of 
the upper level psit) always remains negligible, the pop- 
ulations pi{t) and P2{t) reach the steady state value of 
0.5, and the theoretical curves fit the envelopes of the fast 
oscillating functions. The above results provide a clear 
demonstration of the validity of the two- level approxima- 
tion developed in section II, which is at the basis of the 
subsequent theoretical developments. 
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FIG. 5: Hamiltonian dynamics of the full three-level system 
for the parameters A' = 0.9, g = 0.004, Q, = 0.1, n[ = 0.05, 
Q'2 = 0.1. a) Cavity field mean photon number vs dimension- 
less time: numerical value (solid line) and theoretical value 
(dashed line) . b) Atomic populations: numerical values (solid 
hnes) and theoretical values (dashed lines). 



b) 




5000 At lOOOO 



5000 
At 



10000 



FIG. 7: Full dynamics of the three-level system for the pa- 
rameters as in Fig. O and for k = g^g = 0.00044. a) Cavity 
field mean photon number vs dimensionless time: numeri- 
cal value (solid line) and theoretical value (dashed line), b) 
Atomic populations: numerical values (solid lines) and theo- 
retical values (dashed lines). We used twenty trajectories. 




-5 -5 



FIG. 6: Wigner function for the cavity field state after detec- 
tion of the atom in the ground state 1 1) . The parameters are 
as in Fig. [5] and the dimensionless time is At — 7160. 

VII. CONCLUSIONS 

We have introduced a solvable model of a strongly- 
driven one-atom laser in the optical regime of cavity 
QED. We have shown analytically and numerically that 
the complex dynamics of a three-level atom, dispersively 
coupled to an optical cavity mode and to three laser 
fields, can be well approximated by a two-level atom that 
is resonantly coupled to a cavity mode and a strong co- 
herent field. The effective coupling constant is a com- 
bination of the atom-cavity field coupling constant, the 



amplitude of one of the external lasers and the detun- 
ing parameter. The initial transient regime shows that 
the system is approximately in an entangled atom-cavity 
field state, a Schrodinger cat state, and we show that 
the amount of entanglement depends on the steady state 
mean photon number that is the ratio between the ef- 
fective coupling constant and the cavity decay rate. In 
addition, we propose a scheme for monitoring the whole 
system decoherence based on atomic population measure- 
ments. We find that, for large values of the steady-state 
mean photon number (i.e. in the strong coupling regime), 
and for time larger that the inverse of the effective cou- 
pling constant, the decoherence behavior can be well ap- 
proximated by an exponential decay whose rate is given 
by the effective coupling constant. 

The cavity field subsystem is always in a mixed state 
whose photon number distribution is Poissonian, while 
the atomic subsystem can exhibit coherences. If we mea- 
sure the atomic state at a given time, we can project the 
cavity field in a cat-state with sub-Poissonian or super- 
Poissonian photon statistics depending on the detected 
atomic state. 
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